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Abstract 

We discuss the Yano-Koonin-Podgoretskii (YKP) parametrization of the two- 
particle correlation function for azimuthally symmetric expanding sources. 
We derive model-independent expressions for the YKP fit parameters and 
discuss their physical interpretation. We use them to evaluate the YKP fit 
parameters and their momentum dependence for a simple model for the emis- 
sion function and propose new strategies for extracting the source lifetime. 
Longitudinal expansion of the source can be seen directly in the rapidity de- 
pendence of the Yano-Koonin velocity. 
PACS numbers: 25.75.Gz, 25.75.Ld, 12.38.Mh 
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The two-particle correlation functions C(pi,p2) of identical particles provide direct ac- 
cess to the spatio-temporal evolution of the collision region in heavy-ion collisions. This 
follows from the basic relation [0-|5j between the correlation function C and the emission 
function S(x, K) (here written down for bosons) 

C(K, q) ~l + '^^ A >T |2 . (1) 
V ; \Jd*xS(x,K)\ 2 V ' 

The emission function S describes the phase space (Wigner) density of the boson emitting 
sources, and q = pi — P2, K = (jpi + p2}/2 (with px,p2 being on-shell such that K ■ q = 0) 
correspond to the relative and average 4-momenta of the boson pair. Eq. ([!]) neglects final 
state interactions; for a recent review of methods to include the latter see Ref. ||. All other 
approximations leading to (|l|) are well-controlled [§,0. 

It is the aim of Hanbury-Brown/Twiss (HBT) interferometry to extract information 
about the space-time characteristics of S(x, K) from the measured two-particle momentum 
correlations by "inverting" Eq. (j^). Unfortunately, due to the on-shell conditions for the 
individual particle momenta Pi,pz, this is not possible in a completely model-independent 
way: the time component q° of the relative momentum is constrained by 

n „ „ K 2K K 

,» = /3. q , l3= - = JS - Twr - (2) 



(with E^=\/va 2 + K 2 in the last step @), and thus the inverse Fourier transform of Eq. (U) 
is not unique. 

In practice the analysis of HBT correlation data must therefore be based on a com- 
parison with specific models for the source function S(x,K), with the aim of eliminating 
"unreasonable" models and fitting certain essential parameters (geometric size, freeze-out 
temperature, collective flow velocity, time duration of the particle emission process) in a class 
of "reasonable" model sources. This procedure is enormously simplified by using so-called 
"model-independent" expressions for the HBT parameters ]7|~@] which allow to calculate 
from an arbitrary source function S the characteristic parameters of the two-particle corre- 
lation function C by simple quadrature. 



While it is obvious from its definition that the single-particle spectrum is nothing but 
the zeroth space-time moment of the emission function, 

E *§jC = jd 4 xS{x,K), (3) 

it can also be shown |7|,|9|JTT|] that the two-particle correlation function is essentially deter- 
mined by its (normalized) second order space-time moments. Specifically, to compute the 
correlation function C it is sufficient to approximate the source function S by a Gaussian 
which contains only information on its space-time moments up to second order. To see this 
we write 



--F(K)^(K)F(K) 



SS(x,K), (4) 



S(x, K) = N(K) S(x(K), K) exp 
where 

x»(K) = (x»), x»(K) = x»-x»(K), (B- 1 ) ia ,(K) = (xpX v ), (5) 

with the (K-dependent) expectation values defined as space-time averages over the source 
function 

Then the term SS has vanishing zeroth, first and second order moments and thus contains 
only higher order information on sharp edges, wiggles, secondary peaks, etc. in the source. 
It was shown numerically |TT[ to have negligible influence on the half width of the correlation 
function and to contribute only weak, essentially unmeasurable structures in C(K, q) at large 
values of q. Neglecting SS, the single-particle spectrum @ and the two-particle correlation 
function flU) can be calculated analytically: 

Ek§^ = N(K) det ((iTVCK)) S(S(K), K) , (7) 

C(K,q) = l + exp[-gV(^)( K )] ■ ( 8 ) 

The factor det(B~ l (K)) in (^j can be interpreted JTD| as the generalized 4- volume of the 
emission region for particles with momentum K, V* (K) = det {{x^xj)). However, due 



to the i^-dependent normalization factor N(K), neither V* (K) nor the point a; M (K) of 
maximum emissivity at momentum K can be uniquely unfolded from the single-particle 
spectrum; the latter also drops out from the two-particle correlation function. Only the 
(K-dependent) effective widths ("lengths of homogeneity" |T2|J7[1) (x^x^iK.) of the source 



of particles with momentum K are accessible by HBT interferometry. 

Furthermore, due to the on-shell constraint (^), only 6 linear combinations of the vari- 
ances (xnX u )(K) are actually measurable |fT3[| ; in the case of azimuthal symmetry of the 
source around the beam axis, this number reduces to 4. To make contact with experimental 
correlation data, the redundant components must be eliminated from the exponent of the 
Gaussian in (|8|). It is convenient to do this by using a cartesian coordinate system with z 
along the beam axis and K lying in the x-z-plane. Customarily one labels the z-component 
of a 3- vector by I (for longitudinal) , the x-component by o (for outward) and the y-component 
by s (for sideward). Then from (fj) we see that j3 s = such that 

q° = P±qo + Piqi (9) 

with j3± = \K.±\/K° being (approximately) the velocity of the particle pair transverse to the 
beam direction while fii is its longitudinal component. 

The standard form f?],[|] for the parametrization of the correlation function is obtained 
by using (|9]) to eliminate q° from Eq. (§). One obtains 



C(K, q) = 1 + exp 



i,j=s,o,l 



(10) 



where the 6 HBT radius parameters Rij are defined in terms of the following variances of 
the source function 0,0: 

4(K) = ((x, i,j = s,o,l. (11) 

For an azimuthally symmetric collision region or an azimuthally symmetric sample of col- 
lision events, C(K, q) is symmetric with respect to q s —q s [jnj. Then B 2 ^ = E? sl = 
and 



C(K, q) = 1 + exp [-R 2 s (K)q 2 s - R 2 (K)q 2 - R 2 (K)q? - 2i£(K)g o(ft J , (12) 

with 

R 2 (K) = (y 2 ), (13a) 

R 2 (K) = {{i - (3 J) 2 ) , (13b) 

R 2 {K) = {{z-(3 l i) 2 ), (13c) 
i§(K) = ((5-/3^(5- At)). (13d) 

Clearly these HBT radius parameters mix spatial and temporal information on the source in a 
non-trivial way. Their interpretation in various reference systems, in particular the meaning 
of the generally non- vanishing cross-term R 2 oll was extensively discussed in Refs. [[?|,Pi |lT|p!3 



by analyzing these expressions analytically for a large class of (azimuthally symmetric) model 
source functions and comparing with the numerically calculated correlation function ([!]). An 
important observation resulting from these studies is that the difference 

R 2 iS = Rl - R 2 = Pl(P) - 2(3 ± (xt) + ((x 2 ) - (y 2 )) (14) 

is generally dominated by the first term on the r.h.s. and thus provides access to the lifetime 



At = y (t 2 ) ~ ft) 2 °f the source 0] (more exactly: the duration of the particle emission 
process). However, in heavy-ion collisions, due to rapid expansion of the source one would 
not expect (t 2 ) to be generically much larger than either (x 2 ) or (y 2 ); in the situations 
investigated so far (e.g. it comes out an order of magnitude smaller. In the standard 
fit one is not sensitive to small values of At since Eq. (0) then involves a small difference of 
two large numbers, each associated with standard experimental errors. The factor f3j_ < 1 
in front of (t 2 ) further complicates its extraction, in particular at low K± where At(K) is 



usually largest (see below). Indeed, published experimental results p~5|JT6|1 so far show no 
positive evidence for a finite duration of the particle emission process, in contradiction to 
all physical intuition. 



We will show here that a generalization to azimuthally symmetric systems JT5| , pT| of the 



Yano-Koonin parametrization for a moving source |l8j circumvents this problem [ 19[| . This 



Yano-Koonin-Podgoretskii (YKP) form is based on an elimination in Eq. ([8]) of q a and q s in 



terms of qj_ = Jq 2 + q 2 , q°, and q 3 , using Eq. (Jfty: 
C(K, q) = 1 + exp \-R 2 ± (K) q\ - i?f(K) (q? - (q ) 2 ) - (R 2 (K) + R 2 (K)) (q ■ U(K)f 



(15) 

where U(K) is a (^-dependent) 4-velocity with only a longitudinal spatial component: 

U(K) = 7 (K) (1, 0, 0, i;(K)) , with 7 = • (16) 

V 1 — f 

This parametrization has the advantage that the YKP parameters i?^(K), Rq(K), and 
R 2 (K) extracted from such a fit do not depend on the longitudinal velocity of the observer 
system in which the correlation function is measured; they are invariant under longitudinal 
boosts. Their physical interpretation is easiest in terms of coordinates measured in the frame 
where v (K) vanishes. There they are given by [13| 

Rl(K) = R 2 (K) = (y 2 ) , (17a) 



^(K) = ^-A*j j - ^(f) ~ (z 2 ) , (17b) 
R 2 (K) = ((t- ±x) ) - ^(f) « (?) , (17c) 



where in the last two expressions the approximation consists of dropping terms which were 
found in to be generically small (an extensive and quantitative discussion of this point 
will follow elsewhere 0]). The first expression (|17a ) remains true in any longitudinally 



boosted frame, and we will therefore now concentrate on the other three YKP parameters. 

Eq. (|1 7c|) shows that the YKP parameter Rq(K) essentially measures the time duration 
At(K) during which particles of momentum K are emitted, in the frame were the YKP 
velocity v (K) = 0. The crucial point here is that the smallness of the difference (x 2 — y 2 ) is 
already accounted for directly by the fit, and no potentially small prefactor j3j_ occurs. This 
means that the extraction of Ai(K) from the YKP-parameter Rq(K) is much more direct 
and subject to less statistical uncertainties than in the standard fit. Clearly, this point is 
only true and our suggestive simple spatio-temporal interpretation of the YKP parameters is 



only valid as long as the approximations in (|17D are justified. For realistic emission functions 
they are as we shall show below. 

Eqs. ( |I"7D were written down |13|] in the special frame where v (K) = which we call 
Yano-Koonin (YK) frame. In |13| it was shown that for a large class of models this frame 
essentially coincides with the longitudinal rest frame of the fluid cell around the point x(K) 
of maximum emissivity at momentum K (i.e. the Longitudinal Saddle Point System LSPS 



10|| ) . This was true also for sources which are not longitudinally boost-invariant and for 



which the LSPS and the LCMS (the Longitudinally CoMoving System in which the pion 



pair has (3\ = |14[]) do not coincide. 

We now give model independent expressions, similar to Eqs. (|TT|), (|13|) and ([14]), for the 
YKP fit parameters in an arbitrary observer frame. They are again given in terms of second 
order moments of the source function S(x, K) and thus calculable by simple quadrature. 
The expression for v (K) can then easily be used to establish, analytically and numerically, 
the relationship between the YK frame and the various other frames mentioned above. 

We introduce the following notational shorthands: 

A= ({"-£))• (i8a) 

B=((i-M'), (18b) 



HH) ('-£)>■ 

where £ = x + iy and (y) = (xy) = for azimuthally symmetric sources such that (£ 2 ) = 
(x 2 — y 2 ). In terms of these expressions one finds 



R\ = l - (J (A + B) 2 — AC 2 -A + D ) = D-rC. ( 1 91 > ) 



R 2 = - (A + B) 2 — AC 2 + A — Bj = A-vC, (19c) 
The Yano-Koonin velocity v is zero in the frame where the expression (|18cj) for C vanishes 
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13 ; this fixes also the sign in front of the square root in (|19a|) . For small values of C the 
Yano-Koonin velocity is given approximately by 

C „ {zt) 



A + B (t 2 ) + (P) ' 



(20) 



where in the second approximation we again neglected generically small terms []T3[ propor- 
tional to (zx), (xt), and (x 2 — y 2 ). The accuracy of the approximate expression ( 2Q|) for 



i>(K) was tested numerically and found to be excellent in the situations discussed below. In 
the same limit the expressions for Rq and R 2 simplify to Rq ~ A and R 2 ~ B, in agreement 
with ([Tj 

It is instructive to compare the standard and YKP forms, Eqs. ([12]) and fll5|), for the 
two-particle correlation function. One finds Eq. ( |17a| ) and 

R 2 diS = R 2 -R 2 S = (3 2 ±7 2 (R 2 + v 2 R(j (21a) 
R 2 = (l - A 2 ) R\ + 7 2 (A - vf (Rl + Rl) (21b) 
Rli = P± (-PiRl + 7 2 (A - (Rl + i? 2 )) (21c) 

To invert this set of equations we calculate (cf. Eqs. flI8p) 

A = ^R% s , (22a) 
B = R 2 i ~ -TT-Rli + TprRdm (22b) 



P± ol Pi 

< --^Rli + §R% S - (22c) 
Px Pi 



Inserting this into Eqs. (|II]) gives very cumbersome expressions which provide little physical 
insight. Thus, while the standard HBT radii are easily obtained from the YKP parame- 
ters via ([H]), the converse is not true. This indicates that the YKP parameters are more 
"physical" than the standard HBT radii. Nevertheless, the relations (21) provide a powerful 



consistency check on the experimental fitting procedure of the correlation function, of simi- 
lar value as the relation |i3l , |PT| lim^ ± _^ (i? o (K) — i? s (K)) = which results from azimuthal 



symmetry. 
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We now discuss numerically the dependence of the YKP parameters on the pair mo- 



mentum K. For our study we use the model of Ref. |13[ for a finite expanding thermalized 
source 

S{z,K) = M±C ° SHV ~ Y) (23) 
(27r) 3 ,/27r(Ar) 2 



x exp 



K-u(x) (r-r ) 2 r 2 (v ~ Vof 



T 2(Ar) 2 2R 2 2(Ar?) 2 

Here r = ^x 2 + y 2 , the spacetime rapidity 77 = |ln[(t + z)/(t — z)\ and the longitudinal 



proper time r = \Jt 2 — z 2 parametrize the spacetime coordinates x M , with measure d 4 x 



r dr drjr dr d<p. Y — | ln[(l+/3j) / (1— Pi)] and Mj_ = ym 2 + K\ parametrize the longitudinal 
and transverse components of the pair momentum K. T is the freeze-out temperature, R 
is the transverse geometric (Gaussian) radius of the source, To its average freeze-out proper 
time, At the mean proper time duration of particle emission, and Arj parametrizes 0] the 
finite longitudinal extension of the source. The expansion flow velocity u^{x) is parametrized 
as 

u^(x) = (cosh f] cosh r)t(r), sinh rjt(r) e r , sinh 77 cosher)) , (24) 

with a boost-invariant longitudinal flow rapidity rji = r\ and a linear transverse flow rapidity 
profile 

Vt(r) = Vf (£) ■ (25) 

r]f scales the strength of the transverse flow. Other possible features of the source, like spatial 
and temporal gradients of the freeze-out temperature ||10|| , other freeze-out hypersurfaces or 
different flow profiles, will be discussed elsewhere. 

For the numerical calculations in this letter we have selected one fixed set of source 
parameters: R = 3 fm, to = 3 fm/c, Ar = 1 fm/c, Arj = 1.2, T = 140 MeV. We study only 
pion correlations and set m = m v = 139 MeV/c 2 . Results for different parameter sets as 
well as for kaon correlation functions will be presented in a longer paper |2(| . 
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In Fig. |l]we show the relationship between the YK frame and the LCMS and LSPS. Y is 
the pion pair rapidity (and thus the rapidity of the LCMS), Y YK (Y, K±) the rapidity of the 
Yano-Koonin rest frame, and Y hSPS (Y, K±) the rapidity of the longitudinal rest frame of the 
point x(Y, K±) of maximum emissivity (all rapidities are measured relative to the CMS of 
the source). For pion pairs with large Kj_ both the YK rest frame and the LSPS rapidities 
approach the LCMS rapidity Y, i.e. in this limit all the pions are emitted from a small 
region in the source which moves with the same longitudinal velocity as the pion pair. For 
small K± the YK frame is considerably slower than the LCMS, but faster than the LSPS. 
The linear relationship between the rapidity Y YK of the Yano-Koonin frame and the pion pair 
rapidity Y is a direct reflection of the boost-invariant longitudinal expansion flow. Such a 
behaviour, and thus direct evidence for a strong longitudinal expansion of the source, was 



recently found experimentally in Mg+Mg collisions at 4.4 A GeV/c in Dubna |21] . 

The difference between Y YK and Y LSPS is due to a longitudinal asymmetry of the source 
around the saddle point x(Y, K±); if the source is z-symmetric around x(Y, K±) the YK rest 
frame and the LSPS become identical Both Y YK and Y LgPS exhibit only a very weak 

dependence on the transverse flow of the source; its origin will be discussed quantitatively 

in [gg. 

In Fig. |] we show Rq and R\\ as a function of K± for pion pairs with momentum Y = 



and Y = 3 in the CMS frame and compare these radii with the approximations R\\ ~ J (z 2 ), 



Rq « y (t 2 ) given in Eqs. (|i~7|b,c). The approximation is seen to be exact for vanishing 
transverse flow, rjf = (as already pointed out in ||13]| ). For R\\ it remains rather accurate 
for all K±- values even in the presence of large transverse flow (Fig. |2|d). The parameter 
Ro, on the other hand, is an accurate measure of At only for small K± or sufficiently small 



transverse flow ||13|| . The difference between these two quantities arises from the terms 
-2(xt) /p ± + (x 2 - y 2 ) which were neglected in the second equality of Eq. ( |1 7c| ) . As seen 
in Figs. |2|a,b, these terms can become a serious source of error in the determination of At(K) 
(in our case an overestimate of up to 50% in the most unfavorable case) if the transverse flow 
of the source is very large and not independently known such that it could be corrected for. 
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However, it should be noted |13] that the contamination by these undesired terms is absent 
for pion pairs with small pair momentum K, and that therefore the determination of At 
from the YKP-parameter Rq is particularly clean in the region where its extraction from the 
standard fit according to Eq. (|i4|) is difficult due to the /3^-prefactor. Furthermore, the terms 
that contaminate B% at large K and large transverse flow affect the extraction of At from 
Eq. fll4]) in exactly the same way. This problem can thus not be avoided be selecting either 
the standard or the YKP fitting procedure; by doing and comparing both, in particular 
also for heavier particles, it may be possible to estimate the amount of transverse flow and 
correct for it |20|| . Here it should suffice to say that the associated relative error on At is 
everywhere less than 25% for transverse expansion velocities T]f < 0.3 which we believe to be 
realistic, and that it should decrease for more realistic larger values for the model parameter 
At than the 1 fm/c chosen in Fig. |2]. 



Both the longitudinal region of homogeneity J (z 2 ) and the effective lifetime y (t 2 ) of 
the source decrease for pion pairs with large momenta in the CMS of the source. Asymp- 
totically the effective lifetime becomes equal to the model parameter At = 1 fm/c, but 
low-momentum pions see a much larger value. This is because for low pair momenta the 
longitudinal region of homogeneity R\\ is large, and the correlation function receives also 



contributions from regions freezing out at later times t = ^Tq+z 2 ± At along the surface 
of constant proper time r (z and t measured in the YK rest frame). This is a generic effect 



which should also appear for different source models |22j| . The resulting strong variation of 
(t 2 ) at small K± is again hard to extract from the standard fit because this region is sup- 
pressed by the factor (5\ in Eq. fll4]). Although for large transverse flow Rq does no longer 
exactly trail the lifetime (t 2 ), it clearly reflects this strong K-dependence of the latter at 
small values of K. 

To summarize, we have given model-independent expressions and a detailed physical 
interpretation of the fit parameters for a Yano-Koonin-Podgoretskii fit to the two-particle 
correlation function. We have also established a simple analytical relation between these pa- 
rameters and the "standard" HBT radius parameters which provides a powerful consistency 
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check on the experimental fitting procedure. We clarified the relationship between the YK 
rest frame and the previously introduced LCMS and LSPS frames and argued that the YKP 
fit parameters provide the most intuitive characterization of the local geometric and dynam- 
ical space-time characteristics of the source. An increase of the YK velocity with the pair 
rapidity signals longitudinal expansion of the source. We also pointed out a strong generic 
K-dependence of the effective duration of particle emission which results mainly from the 
fast longitudinal expansion of the source, but is also modulated by transverse expansion. 
We hope that all these predictions will soon be checked experimentally in relativistic heavy 
ion collisions. 

This work was supported by grants from DA AD, DFG, NSFC, BMBF and GSI. We grate- 
fully acknowledge discussions with H. Appelshauser, S. Chapman, D. Ferenc, M. Gazdzicki, 
and P. Seyboth. 
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FIGURES 

FIG. 1. (a) The rapidity of the YK frame as a function of the pion pair rapidity Y (both 
measured in the CMS frame of the source), for various values of the transverse momentum K± 
of the pair and two values for the transverse flow rapidity rjf. (b) Same as (a), but shown as a 
function of K± for different values of Y. The curves for negative Y are obtained by reflection 
along the abscissa, (c) The difference Y YK — Y LSPS between the rapidity of the YK frame and the 
longitudinal rest system of the saddle point, plotted in the same way as (a), (d) Same as (c), but 
shown as a function of K± for different values of Y . 

FIG. 2. (a) Rq and yj (t 2 ) as a function of M± for three values of the transverse flow rapidity 
rjf, for pion pairs with rapidity Y = in the source CMS frame. The lifetime \J (t 2 ) is evaluated 
in the YK rest frame (which in this case coincides with the CMS frame), (b) Same as (a), but for 
pions with rapidity Y = 3 in the CMS frame, (c) and (d): Same as (a) and (b), but for i?y and the 
longitudinal length of homogeneity y/ (z 2 ) in the YK rest frame. For Y = 0, i?y and \J (z 2 ) agree 
exactly because = in the YK frame. 
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